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Abstract. - We study the forced fluid invasion of an air-filled model porous medium at 
constant flow rate, in f +f dimensions, both experimentally and theoretically. We focus on the 
non-local character of the interface dynamics, due to liquid conservation, and its effect on the 
scaling properties of the interface upon roughening. Specifically, we study the limit of large 
flow rates and weak capillary forces. Our theory predicts a roughening behaviour characterized 
at short times by a growth exponent I3\ = 5/6, a roughness exponent ai — 5/2, and a dynamic 
exponent zi = 3, and by /32 — 1/2, 02 — 1/2, and Z2 — 1 ai long times, before saturation. This 
theoretical prediction is in good agreement with the experiments at long times. The ensemble of 
experiments, theory, and simulations provides evidence for a new universality class of interface 
roughening in 1 -I- 1 dimensions. 



Studies on the morphology of interfaces moving in disordered media under nonequilibrium 
conditions constitute an active field of research |Q] . One relevant example is the forced invasion 
of an air- filled porous medium by a viscous wetting fluid such as oil or water, which gives rise 
to a nonequilibrium rough interface . The roughening process of an initially flat interface 
is described in terms of the intcrfacial root- mean-square width w. In many systems w follows a 
Family- Vicsek dynamical scaling |7j: w ^ t'^ ioi- t < t^ , w ^ ioi- t > t^ and <x ^ L^, with 
a = zp. Here t^ is a saturation time, (3 the growth exponent, a the roughness exponent, and z 
the dynamic exponent. The roughness exponent can also be obtained from the power spectrum 
S{q,t), which is less sensitive to finite-size effects. We have g S{q,t) = q-^^°'+'^lsFviqt^^''), 
where spv obeys spviu) constant when u ^ 1, and spviu) ~ m^""*"^ when u <^ 1. 

Experiments of forced fluid invasion (FFI) in air-fllled packings of glass beads, in 1+1 
dimensions, have given roughness exponents in the range a — 0.6 — 0.9 pHq|. The dispersion 
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Fig. 1 - Several consecutive close-up views of a piece of the oil-air interface, taken at 150 s time 
intervals. The flow rate is Q = 10 ml/h. 

reflects that the effective values of a, determined by measuring widths over several scales 
at saturation, depend on the capillary number Ca 0- On the other hand, there are few 
experimental determinations of the exponent /3, because the large relative strength of capillary 
to viscous forces made the growth regime extremely short in most experiments. 

In this letter we report new experimental results of dynamic interfacial roughening in FFI, 
in 1-1-1 dimensions, at constant injection rate. In our case the model porous medium is a Hele- 
Shaw cell with controlled spatial fluctuations of the gap thickness. While the viscous pressure 
field and the surface tension in the plane of the cell keep the interface smooth on large length 
scales, fluctuations in the gap thickness produce local fluctuations in capillary pressure which 
roughen the interface. Thus the governing physics is the same than in a real two-dimensional 
porous medium, and we have the possibility of controlling the statistical properties and the 
relative strength of the disorder. 

In our experiments a silicone oil |^ displaces air in a horizontal Hele-Shaw cell, 190 x 550 
mm^, made of two glass plates 19 mm thick. Fluctuations in the gap thickness are provided 
by a fiber glass substrate which contains a random distribution of square copper islands filling 
35% of the substrate. Each island is 0.06 ± 0.01 mm thick and has 1.500 ± 0.005 mm lateral 
size (Fig. |]). The gap spacing from the substrate to the top plate is 5 = 0.36 ± 0.05 mm. 
Since the maximum width of the meniscus is only one half the gap width, the interface can 
be considered one-dimensional at the length scale of the copper islands. The oil is injected 
at one side of the cell using a constant flow-rate syringe pump. A flat front is first prepared 
by keeping the oil at rest in a transverse copper track separated 2 mm from the beginning of 
the noise. Initially the oil is pushed at high flow rate (300 ml/h) for about 4 s to keep the 
interface smooth while it enters the noise. Next, at a time defined as i = 0, the fiow rate is 
set to its nominal value Q, and the oil- air interface is monitored as a function of time using 
two CCD cameras, with a final spatial resolution of 0.37 mm/pixel. This imposes a cut-off 
at small length scales, about twice the maximum width of the meniscus. We choose a range 
of fiow rates (10 < Q < 100 ml/h) such that the contact line overcomes the copper islands 
rather easily, without pinning or forming overhangs. For Q = 10 ml/h (Fig. |l|) the average 
interface velocity is w = 0.038 ± 0.005 mm/s, which we take as a reference. 

Our results for w(t) are presented in Fig. ||(a). As expected, both the saturation time and 
the interfacial width w at saturation depend on the average velocity at which the interface 
is driven. We observe also large fluctuations in the interfacial width, consistent with the 
observations reported in Q]. The fluctuations in w are highly reproducible in experimental 
runs with the same realization of the disorder. Disregarding fluctuations, the growth of the 




Fig. 2 - a) Double-logarithmic plot of w{t) for four different interface velocities. Each curve is an 
average over 6 experiments (3 runs on 2 different realizations of the noise). In each experiment ■w'^{0) 
has been substracted from iiP{t), in order to minimize the effect of the initial condition. The error 
bars give a measure of the overall dispersion between different experiments, b) Power spectrum of 
the interface profiles at successive time intervals, for the experiments at velocity v. The vertical line 
indicates the spatial scale of the disorder. 



interfacial width is consistent with a power-law w ^ with (3 — 0.47 ±0.04, extending about 
three orders of magnitude in t (for measurements at velocity v) until the width saturates. 

The power spectrum for the set of measurements at velocity v is shown in Fig. ^(b). At 
short times the spectrum displays a plateau for small q and a power law decay at large q 
as expected. As the interface advances in the disorder the value of the plateau increases, 
and we observe the emergence of another scale-independent (power-law) behaviour for an 
increasingly growing range of spatial scales, in agreement with the scaling of S(q, t) expected 
for u = qi^/^ ^ 1. This behaviour is also observed at velocity 2v. The measured exponents are 
-2.2±0.2 (a = 0.6±0.1) at small q (long length scales) and -2.8±0.2 (a = 0.9±0.1) at large q 
(short length scales) . The spectra of the experimental interfaces have been calculated imposing 
periodic boundary conditions. To this purpose, we have substracted from the interface the 
straight line joining its two ends. This correction eliminates an artificial slope that is specially 
important at saturation, where the difference in height between the two ends is maximum pO| . 

One interesting aspect of roughening in fluid flows is the nonlocal nature of the interfacial 
dynamics, due to liquid conservation. The importance of nonlocality in this problem has 
been already pointed out by a number of authors |4[ ^, |ll|] . Very recently this issue has been 
addressed explicitly in the theoretical formulation of FFI at constant pressure [|l2| and of 
spontaneous imbibition In this letter we address the same issue for FFI at constant 

flow rate. This driving condition has not been explored theoretically yet, in spite of being 
often used in experiments. 

Our model is based on a time dependent Ginzburg Landau model with conserved order 
parameter (model B in Ref. [Q): 

^ = VMV(-0 + 03_,2v2^). (1) 



(j) is the order parameter with equilibrium values (j)eq = ilj and e is the interfacial width. M 
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is the mobility, which we consider fluctuating in space: 

M(^,r) = |^(l+^W) J--^ (/>2 (2) 
^ [0 for air (0 < 0) ^ ' 

with X a macroscopic mobility and ^(r) a weak static disorder with a spatial correlation: 

< e(r)C(r') >= 2DC (^^) , (3) 

where A is the disorder correlation lenght and C is normalized. Dube et al. have proposed 
a similar formulation for spontaneous imbibition. In their problem the interface is driven only 
by capillary forces, and moves with average velocity v '-^ Here, instead, nonlocal viscous 

forces are dominant over the local fluctuations of capillary pressure: our interface is driven 
by a constant flux 70, which gives rise to a steady state rough interface moving with constant 
noise-dependent v. The disorder is also introduced differently in the two approaches. In our 
two-dimensional model the fluctuations in gap thickness are represented (phenomenologically) 
by a fluctuating mobility. 

The macroscopic limit of the model is obtained by an asymptotic expansion in orders of 
e, following a procedure described in Ref. p^ : 

V(l+e(r))VP = , 

P = -Tk , (4) 

where the pressure field in the bulk of the liquid is given by P = {(j){r,t) — (j)f.q)/{2(j)eq), the 
surface tension is given by F = (20eg)~^ / dy{dy(j)stY , 4>st is the steady planar solution of Eq. 
(|l|), K is the curvature of the interface, and w„ its normal velocity. £,{x,yint) is taken at the 
interface. The deterministic part of Eq. reproduces the macroscopic equations for oil-air 
displacements in a Hele-Shaw cell [ p^ . 

The interfacial equation in Fourier space can be derived from Eq. using Green function 
analysis. We define P = P + jo/K{y — jot), where 70 = 7o/(40eg). In this case, the Green 
function G{r,r') obeys V^G — 6{r — r'). The interfacial equation for the local deviations of 
the interface from its mean position, h{x, t) = yint — vt, reads (to first order in ^): 



ix'"n{s') ~ lofiy _ 1 



K 2 



K 



+ / ds'hVGis,s') -rK{s') + ^{h{s') + {v-^o)t) + / ds'G{s,s')^as\y^nt), (5) 



K 



K 



where s is the contour variable on the interface and h and y are the unitary vectors perpendic- 
ular to the interface and along the y direction respectively. We have neglected volume terms 
which do not contribute to the scaling. A distinguishing feature of Eq. (||) is the presence 
of integral terms, which account for the nonlocal character of the interfacial dynamics. The 
deterministic part is equivalent to the equation for unstable displacements (air displacing oil, 
70 < 0) derived in Refs. 0,11). 

In this problem, due to the presence of a correlated disorder, two different time scales of 
dynamic origin are present. One of them is related to the ballistic dynamics of the interface 
due to the driving and corresponds to the time in which the interface advances the noise 
characteristic length in the y-direction, to = A/70. The other temporal scale corresponds to 
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the dynamics that couples different points of the interface. For scales not much greater than A 
this dynamics is diffusive. The associated temporal scale is the time required by the diffusion 
along the interface to reach distances of the scale of the disorder, and is given by tu = }? /K . 
Here, we consider to ^ tn- This condition is fulfilled if A ^ K/'Jq. In this limit, the noise 
appears as persistent in the y-direction at scales comparable to A. Therefore, we will take it 
as equivalent to a columnar (only x-dependent) noise. In this case, the stationary velocity 
= 7o + {iq/L) J £,{x)dx, where L is the system size in the x direction. 
In Fourier space, for small deviations from a flat interface {\q\h <^ 1) we have: 



dhq 

dt 



(6) 



We observe that short wavelength fluctuations are damped by surface tension (first term of 
the r.h.s.), while long wavelength fluctuations are damped by the advancement of the front, 
driven by the external flux (second term of the r.h.s.). Eq. (|6|) leads to two different temporal 
regimes. The early growth regime is dominated by the first term of (|^), and the dynamic 
exponent is zi = 3. By linear scaling analysis we obtain z — a — 1/2 = 0, which gives ai = 5/2 
and Pi = 5/6. Since ai > 1, the interface is superrough in the early time stages. At longer 
times the second term of (^ dominates the dynamics, and we get Z2 — 1, a2 — 1/2, and 
P2 — 1/2. The characteristic crossover time between the two growth regimes is given by 
^ Cc' ~ (KT/^of'. 
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Fig. 3 - a) Temporal evolution of the interfacial width, from the numerical integration of Eqs. (|l])-(P), 
for a system of size L = 256. b) Power spectrum at different time intervals, for the same case of (a). 
The vertical line is the spatial scale of the disorder. The interfacial width at saturation vs. the system 
size L is shown in the inset. 



We have checked these predictions by numerical integration of Eqs. (|i|)-(§) in a rectangular 
lattice with periodic boundary conditions in the x-direction. A constant flux 70 has been 
maintained by flxing the value of (/i at a certain distance behind the interface (typically 20 
space units). To mimic the experiments, we have assumed that the oil mobility can take two 
different values, so that ^ is a dichotomous noise of values zLD. The noise is defined in boxes 
of side Id = A space units, randomly distributed with probability 0.35 for the positive value. 
The values of Id and 70 satisfy the requirement of noise persistence in the y direction. Our 
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results for 70 = 0.05, K = 1, e — 1 and D — 0.5, averaged over 25 realizations of the noise, 
are presented in Figs, ^(a) and ||(b). 

Fig. |^(a) is a log- log plot of u'(t), including the analytical prediction for the growth expo- 
nents. The numerical results are consistent with a crossover between /3i = 5/6 at short times 
and (32 — 1/2 at longer times. Fig. |3|(b) shows a log- log plot of S{q,t). At short times the 
spectrum displays a plateau for small q, followed by a power law corresponding to ai ~ 5/2 
for large q. At longer times, S{q, t) shows a data collapse at small q into a power law corre- 
sponding to a2 — 1/2. The same results are obtained from the analysis of w vs. system size 
L, shown in the inset of Fig. |^(b). We observe the predicted crossover between the two growth 
regimes. Using the scaling relation a = zf3, the dynamic exponents are zi ~ 3 at short times 
and Z2 ~ 1 at long times, in agreement with the analytical results. 

Although we have experimental indications of very fast growth in the earliest time stages, 
we have not been able to identify experimentally the short-time scaling regime predicted by 
the model. This regime is difficult to access experimentally because it is very short, it is 
obscured by the transient originated by setting Q to its nominal value, and, in addition, the 
short-time behaviour of w{t) is strongly dependent on the definition oi t = 0. Nevertheless, 
the experimental power spectrum shows that the short length scales saturate with a larger 
roughness exponent than the long length scales, which points to a different mechanism at the 
two scales, in the same direction as the model. The experimental and calculated values of the 
roughness exponent at short length scales do not coincide because the details of the physical 
mechanisms operative at the shortest length scales (capillary phenomena) are not properly 
captured by the model. 

Concerning the long time and long length scaling regime, several authors have argued that 
FFI should fall in the KPZ universality class {(3 = 1/3 and a = 1/2 in 1-1-1 dimensions) in the 
limit of large Ca considered here |9| . Our experimental and theoretical results contradict 
this prediction for the exponent (3, which is consistently found to be /3 ~ 1/2. Concerning 
the roughening exponent a, however, previous experimental measurements ^ and our own 
experiments give a > 0.6 or larger. Although the value of a decreases with increasing Ca, the 
theoretical limit a = 0.5 is presumably not accessible because the saturation width w falls to 
magnitudes comparable to the gap thickness or to the typical pore size as Ca becomes very 
large. 

In conclusion, we have studied the dynamics of interfacial roughening in forced fluid inva- 
sion of a model porous medium at constant flow rate, in 1-1-1 dimensions. We have focused on 
displacements in which the viscous pressure field dominates over the fluctuations in capillary 
pressure (large Ca, weak disorder), and found a new universality class with two distinct time 
regimes. 
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